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ABSTRACT 

We present a comprehensive spectral analysis of all INTEGRAL data obtained so far for the 
X-ray-bright Seyfert galaxy NGC 4151. We also use all contemporaneous data from RXTE, 
XMM-Newton, Swift and Suzaku. We find a linear correlation between the medium and hard- 
energy X-ray fluxes measured by INTEGRAL, which indicates an almost constant spectral 
index over six years. The majority of INTEGRAL observations were made when the source 
was either at a very bright or very dim hard-X-ray state. We find that thermal Comptonization 
models applied to the bright state yields the plasma temperature of =* 50-70 keV and its op- 
tical depth of =b 1.3-2.6, depending on the assumed source geometry. For the dim state, these 
parameters are in the ranges of =* 180-230 keV and =* 0.3-0.7, respectively. The Compton 
parameter is y a 1 for all the spectra, indicating a stable geometry. Using this result, we can 
determine the reflection effective solid angles associated with the close and distant reprocess- 
ing media as =* 0.3 x 2n and 0.2 x 2n, respectively. The plasma energy balance, the weak 
disc reflection and a comparison of the UV fluxes illuminating the plasma to the observed 
ones are all consistent with an inner hot accretion surrounded by an outer cold disc. The disc 
truncation radius can be determined from an approximate equipartition between the observed 
UV and X-ray emission, and from the fitted disc blackbody model, as ~ 15 gravitational radii. 
Alternatively, our results can be explained by a mildly relativistic coronal outflow. 

Key words: accretion, accretion disks - galaxies: active - galaxies: individual: NGC 4151 - 
galaxies: Seyfert - gamma-rays: observations - X-rays: galaxies 



1 INTRODUCTION 

Continuum properties of the hard X-ray and soft y-ray emission 
from radio-quiet Seyfert galaxies nuclei are relatively well known 
thanks to the many satellites operating during the last decades. 
The spectra are commonly approximated by a phenomenological 
e-folded power-law model, F(E) oc E l ~ r exp(—E/E c ) (where T is 
the photon index and E c is the e-folding, or cut-off, energy), ac- 
companied by a Compton reflection component. A recent study 
based on a local (z < 0. 1) sample of 105 Seyfert galaxies observed 
with BeppoSAX pr esents the aver age parameters of this model for 
Seyfert Is and 2s JDadinj[20o3) . The mean T is a 1.9 (1.8) for 
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Seyfert 1 (2) nuclei, E^ is a 230 (380) keV, and the average relative 
strength of reflection, R = Q/2n (where Q. is the effect ive solid an- 
te su btended by the reflector), is a 1.2 (0.9). See lBeckmann et al.l 
2009) for the similar INTEGRAL results. 



This phenomenological model is too crude for studying of 
the physics of the central engine of Seyferts. There is a general 
consensus that their Xy photons are from (predominantly ther- 
mal) Comptonization by hot electrons of some seed soft photons. 
The seed photons may come from an optically-thick accretion disc 
or clouds in the vicinity of the hot plasma, or may be inte rnally 
produced by the synchrotron process, e.g., Ixie et al.l ( 120101) . For 
example, non-simultaneous data for the Seyfert NGC 5548 from 
ROSAT, Ginga and CGRO/OSSE have yielded the plasma tem- 
perature, kT e a 55 ke V, and the Thomson optical depth, r a 2 
The average OSSE spectra of 11 Sy Is 
70-80 keV, 



( Mag dziarz et alj 11998 

and 8 S y 2s have yielded similar parameters, kT e 

T — 1.7 dZdziarski. Poutanen. & JohnsonluOOfJ) . On the other hand, 
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IPetrucci et al.l ( 120001) obtained a much higher kT s m 250 keV (and 
t =2 0.2-0.4) for NGC 5548 using BeppoSAX data and a differ- 
ent anisotropic Comptonization model in the slab geometry. Then, 
IPetrucci et al.l ( 12001 at) , hereafter P01, found, with the same model, 
a similarly high values of kT e =* 170-320 keV (and t = 0.05-0.20) 
for 6 Sy Is observed by BeppoSAX. We note, however, that such 
high kT s yield spectra being well above the OSSE fluxes at > 100 
keV for either NGC 5548, NGC 4151 or the average Seyferts. 

The nearby, z = 0.0033, Seyfert 1.5 galaxy NGC 4151 is the 
second (after Cen A) brightest persistent AGN in the 20-100 keV 
band. Its electromagnetic spectrum has been extensively studied, 
with its properties in the radio, infrared, optical, ultraviolet and 
Xy bands well established. Its black-hole mass h as been estimate d 
based on reverberation as (4.6+°*) x 10 7 M© teentz et al1l2006h . 
Although this estimate is subject to a systematic uncertainty by a 
factor of ~3-4, it is consistent wit h the estimate based on scaling of 
the power spectrum to Cyg X-l o flCzernv et alj ( 120011) . The overall 
picture of the region surrounding the nucleus, based on HST/STIS, 
Chandra and XMM-Newton data, is that the complex absorption 
observed for NGC 4151 is due to massive outflows, presumably 
disc winds, forming several regions charact erized by a range of 
column densities, A^ h, and ionization levels dSchurch et alj|2004c 
iKraemer et alj 120061) . This explains earlier results from ROSAT, 
ASCA, Ginga and RXTE satellites, with the medium-energy X- 
ray spectra fitted only after applying a complex absorber model, 
consisting of several c omponents fully and partially covering 
the central source (e .g., IZdziarski. Johnson. & Magdziarzlll996al : 
IZdziarski et alj2002l hereafter Z02). The a bsorber undergoes rapid 
daily changes of ,/Vh dPuccetti et alj |2007). Also, a narrow Fe K a 
line with the equivalent width of 50-2 00 eV is observed (Z02; 
ISchurch et alj2003l ; |Pe Rosa et alj2007l) . 

NGC 4151 has frequently been observed by CGRO and 
Bepp oSAX. The e-folded power law model y ields E c = 50-200 
ke V d Johnson etalll 19971: |Pe Rosa et al]|2007 ). Thermal Compton 



models yield kT t a 50-80 keV, t a 1-2 djohnson et alj|l997l : 
Z02) . The strength of Co mpton reflection is moderate, R ss 0.4 
(e.g jDe Rosa et alj2007l) . 

NGC 4151 was observed by INTEGRAL 8 times as a primary 
or secondary target. The r esults of the firs t dedic ated observation 
of 2003 May are given in lBeckmann et alj 1(2005). The object was 
found in a spectral state similar to those observed with OSSE, but at 
the highest flux level ever noticed. There was no significant spectral 
variability despite ~ 50 per cent variations of the flux. A thermal 
Compton model fitted to the summed spectrum yielded kT c a 90 
keV, t a 1.3, R a 0.7, JV H -7x 10 22 cm" 2 . Puring almost all later 
INTEGRAL observations the source was found at moderate to very 
low flux levels. This has allowed us to study the spectral proper- 
ties of the dim state of this AGN with unprecedented precision and 
compare them with the bright state of 2003 May. 



2 OBSERVATIONS AND DATA REDUCTION 

We use all NGC 4151 data collected by INTEGRAL as of 2010 Jan- 
uary, see TableQ] Good quality data come from the dedicated NGC 
4151 observations made in 2003 May, 2007 January, May and De- 
cember, and 2008 May. The remaining data are taken from observa- 
tions of the Coma cluster, Mrk 273, NGC 4736, M5 1 and Mrk 42 1 , 
when NGC 4151 was almost always seen at an off-axis angle of 
9°-15°. We select the data with the off-axis angle < 15° for ISGRI 
and SPI, and < 3° for the JEM-X and OMC. The data have been 
reduced using the Offline Scientific Analysis (osa) 7.0 provided by 



the INTEGRAL Science Data Centre dCourvoisier et alj2003l) . with 
the pipeline parameters set to the default values. Since the osa v. 8.0 
of 2009 August brought a major improvement for the JEM-X data 
analysis, we have used it for these data. The ISGRI and SPI spectra 
and light curves have been extracted with the standard spectral ex- 
traction software, including the catalogue sources NGC 4151, NGC 
4051, NGC 4138, Mrk 766, NGC 4258, NGC 4395, Coma Clus- 
ter, NGC 5033, Mrk 421 and Mrk 268. The JEM-X spectra have 
been obtained from mosaic images stacked for a given observation. 
We use the standard response files for all instruments. T he OMC 
magnitude was converted into flux using the calibration of ljohnsonl 
dl966l) , F(V=0) = 3.92 x 10-' erg cm 2 s~' A" 1 . 

To better constrain the spectra at low energies, we supple- 
mented the INTEGRAL data by all available X-ray observations of 
NGC 4151 taken since 2003 by RXTE, XMM-Newton and Suzaku, 
see Table|2] The RXTE PCU2 light curves and spectra and HEXTE 
spectra were extracted with the heasoft 6.5.1, using standard selec- 
tion criteria. The 5 XMM-Newton observations (denoted here XI- 
X5) were analyzed with sas 8.0.1. We used only the EPIC pn data, 
and excluded periods of high or unstable background. The spectra 
XI and X2 have been found to be compatible with each other; we 
have therefore added them together. For the single Suzaku obser- 
vation, the data in the 3x3 mode were reduced with heasoft 6.6. 
The spectra of the front-illuminated CCD, XIS0 and XIS3, were 
added together, whereas the spectrum of the back-illuminated de- 
tector XIS 1 was used separately. We also extracted the HXD/PIN 
spectrum using the standard procedure. Although there have been 
a number of Chandra NGC 4151 observations, all those without 
gratings that were public at the time we started the analysis were 
from before 2000 March. 

Tables [T}j2] identify the data used for our three main spectral 
sets, bright (B), medium (M), and dim (D), whose selection is based 
on ISGRI flux. The remaining INTEGRAL data, not marked with 
a letter in Table \T\ were not used for spectral analysis. We have 
excluded them because either there were no corresponding high- 
quality spectra from X-ray satellites (e.g., Revs. 0809-081 1) and/or 
their flux level was outside the range assumed by us (e.g., beginning 
of Rev. 0074 and Rev. 0076 in the case of the B state). The data 
sets from the other satellites were not always simultaneous with 
the INTEGRAL data. Thus, we assign them to one of the three sets 
based again on the flux, using the monitoring data from Swift/BAT 
(see Sections|3l l4.1| for details of the selection). 



3 VARIABILITY 

NGC 4151 has been observed with Xy detectors for almost 40 
ye ars. The satel l ite an d balloon results up to 1988 were compiled 
bv lPerotti et alj dl99ll) . Then, it has been observed above 20 keV 
by Ginga, GRANAT, CGRO, BeppoSAX, INTEGRAL, Swift and 
Suzaku. Fig. [TJa) shows the 20-100 keV flux observed since 1972 
October compared to the optical flux and the RXTE/ASM count 
rate. Because the spectra from earlier observations are not publicly 
available, the 20-100 keV fluxes were determined using the values 
of the flux at 35 keV and of the photon index presented in table 2 
of IPerotti et alj dl991l) . For this reason, we do not show their un- 
certainties . Similarly, the GRANA T fluxes were determined using 
table 3a o flFinoguenov et alj dl995ri . where the ART-P and SIGMA 
spectra are fitted by a power-law model. For all other observations, 
we use the spectra from heasarc (CGRO/OSSE, BeppoSAXfPDS) 
or spectra extracted by us (INTEGRAL/lSGKl, SuzakufPIN), where 
the 20-100 keV flux was computed from a power-law model fit 
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Table 1, The observation log of INTEGRAL. The ISGRI exposure is effective, corresponding to fully coded observations, and the JEM-X one is summed over 
all pointings with the off-axis angle < 3°. The last column gives the spectral set as defined in Section FTTl B = bright, D = dim, M = medium. 



Revolution 


Start time, UTC (MJD) 


End time, UTC (MJD) 


Eff. exposure [s] 


JEM-X exposure [s] 


Spectral set 


0036 


2003-01-29 16:43 (52668.697) 


2003-01-31 08:18 (52670.346) 


11296 


- 


- 


0071 


2003-05-14 12:42 (52773.529) 


2003-05-16 11:39 (52775.485) 


9908 


- 


- 


0072 


2003-05-17 13:22 (52776.557) 


2003-05-18 19:56(52777.831) 


5376 


- 


- 


0073 


2003-05-20 18:44 (52779.781) 


2003-05-22 10:23 (52781.433) 


20441 


- 


- 


0074 


2003-05-23 08:35 (52782.358) 


2003-05-25 11:42 (52784.488) 


117528 


91633 


B 


0075 


2003-05-25 21:03 (52784.877) 


2003-05-28 05:33 (52787.231) 


153481 


153481 


B 


0076 


2003-05-28 20:51 (52787.869) 


2003-05-29 12:13 (52788.509) 


39533 


39533 


- 


0274 


2005-01-10 16:05 (53380.670) 


2005-01-12 17:36 (53382.733) 


9146 


- 


- 


0275 


2005-01-13 19:36 (53383.817) 


2005-01-15 04:20 (53385.181) 


11745 


- 


- 


0310 


2005-04-29 16:32 (53489.689) 


2005-04-30 03:46 (53490.157) 


1734 


- 


D 


0311 


2005-05-01 08:54 (53491.371) 


2005-05-02 20:28 (53492.853) 


3290 


- 


D 


0312 


2005-05-03 22:44 (53493.947) 


2005-05-04 09:57 (53494.416) 


1062 


- 


D 


0317 


2005-05-19 18:04 (53509.753) 


2005-05-21 01:09 (53511.048) 


10144 


- 


D 


0318 


2005-05-22 01:03 (53512.044) 


2005-05-24 06:37 (53514.276) 


7894 


- 


D 


0324 


2005-06-09 17:55 (53530.747) 


2005-06-10 23:12 (53531.967) 


9748 


- 


D 


0448 


2006-06-14 09:54 (53900.413) 


2006-06-17 00:03 (53903.002) 


29525 


- 


D 


0449 


2006-06-17 10:19 (53903.430) 


2006-06-19 23:47 (53905.991) 


28497 


- 


D 


0450 


2006-06-20 09:23 (53906.391) 


2006-06-22 23:29 (53908.978) 


26974 


- 


D 


0451 


2006-06-23 09:06 (53909.379) 


2006-06-25 23:10 (53911.965) 


26844 


- 


D 


0521 


2007-01-18 16:45 (54118.698) 


2007-01-20 23:28 (54120.978) 


122878 


122878 


D 


0522 


2007-01-21 16:35 (54121.691) 


2007-01-24 06:07 (54124.255) 


123751 


105848 


D 


0561 


2007-05-18 09:31 (54238.397) 


2007-05-20 16:11 (54240.674) 


110696 


1 10696 


D 


0562 


2007-05-22 00:44 (54242.031) 


2007-05-23 21:34 (54243.899) 


100729 


78177 


D 


0563 


2007-05-24 09:06 (54244.379) 


2007-05-25 06:31 (54245.272) 


46903 


46903 


D 


0634 


2007-12-22 17:53 (54456.745) 


2007-12-25 07:19 (54459.305) 


131811 


131811 


M 


0636 


2007-12-28 17:26 (54462.726) 


2007-12-31 01:53 (54465.078) 


119738 


100028 


M 


0678 


2008-05-02 21:17 (54588.887) 


2008-05-04 22:36 (54590.942) 


24805 


- 


M 


0679 


2008-05-05 20:47 (54591.866) 


2008-05-07 07:09 (54593.298) 


17266 


- 


M 


0809 


2009-05-31 11:31 (54982.480) 


2009-05-31 20:58 (54982.874) 


20019 


2891 


- 


0810 


2009-06-01 05:49 (54983.242) 


2009-06-03 16:10(54985.674) 


112086 


46546 


- 


0811 


2009-06-04 05:37 (54986.234) 


2009-06-06 10:24 (54988.433) 


108501 


40239 


- 



Table 2. The observation log for the other X-ray satellites. The exposure times are given for the RXTE PCU2 (where the number of pointings added together 
is shown in parentheses), EPIC pn (XMM-Newton) and XIS (Suzaku) detectors. The last column shows the data set as defined in Table[T] and, in the case of 
XMM-Newton, its consecutive observation number. 



Obs. ID 



Start time, UTC (MJD) 



End time, UTC (MJD) 



Exposure Spectral set 



80416-01-01-(00-09) 

92113-08-(06-09) 

92113-08-(22-35) 



2003-05-24 05:49 (52783.242) 
2006-05-13 02:28 (53868.103) 
2006-12-23 12:51 (54092.535) 



RXTE 
2003-05-29 04:18 (52788.179) 
2006-06-23 15:42 (53909.654) 
2007-06-22 08:57 (54273.373) 



15760(11) 


B 


4240 (4) 


D 


14272 (14) 


D 


5821 


XI, B 


11389 


X2, B 


12291 


X3,B 


27980 


X4,D 


21468 


X5,M 



0143500101 
0143500201 
0143500301 
0402660101 
0402660201 



XMM-Newton 



2003-05-25 01:38 (52784.068) 
2003-05-26 20:35 (52785.858) 
2003-05-27 15:17 (52786.637) 
2006-05-16 06:22 (53871.265) 
2006-11-29 17:20 (54068.722) 



2003-05-25 06:54 (52784.288) 
2003-05-27 01:51 (52786.077) 
2003-05-27 20:33 (52786.856) 
2006-05-16 17:35 (53871.733) 
2006-11-30 08:00 (54069.333) 



701034010 



2006-12-18 20:05 (54087.837) 



Suzaku 
2006-12-21 09:14 (54090.385) 



124980 



(which was done in the 50-150 keV range for the OSSE data) 
and the uncertainty was determined from the relative error of the 
summed count rate in the fitted band. Fig.[TJb) shows the 20-100 
keV fluxes from the CGRO/BKISE in 7-d bins (together with those 
from OSSE ). They have been obtained from the 20-70 keV fluxes 
given by IParsons et al.l d 199a) by multiplying by 1.1 (assuming 



T = 1.8) and then dividing by the nor malization factor with respect 
to OSSE of 1.42 JParsonset alii 19981) . 

Fig. EM gives the 1.5-12 keV count rate of RXTE/AS^ 
in 30-d bins. The blue curve in Fig. [TJc) shows the opti cal data 
from the Crimean observatory in the V band (5500 A. ICzernv et al.l 



http://xte.mit.edu/asmlc/ASM.html 
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g I 10 




2003 2004 2005 2006 2007 2008 2009 



43000 44000 46000 48000 50000 52000 54000 
Time [MJD] 

Figure 1. Light curves ofNGC4151 since 1972. (a) The 20-100 keV flux 
from the compilation by Perotti et al. 1 1991), magenta diamonds, and from 
GRANAT/SIGMA (cyan squares), CGRO/OSSE (red filled circles), Beppo- 
SAX/PDS (green triangles), INTEGRAL/ISGR1 (blue dots), and Suzaku/PIN 
(the orange open circle). The dotted lines show the extrema of the 20-100 
keV flux from ISGRI. (b) The 20-100 keV flux from CGRO/BATSE (green 
error bars) together with that from OSSE (red filled circles; sa me as in a), (c) 
The R XTE/ ASM count rate in 30-d bins, (d) The 5500 A flux iCzernv et alj 
120031 blue lines), and from I NTEGRAL/OMC (vertical blue bars at MJD 
> 52000), and the 5117 A flux IShapovalova et al.l2 008. black dots). 



2003), which, after MJD 52000, are supplemented by the INTE- 
GRAL/OMC data. Since the V band data do not overlap much in 
time with the ASM data, w e also show the 5 1 17 A flu xes from var- 
ious observatories given bv lShapovalova et alj (12008). 

The 1.5-12 keV count rate appears well correlated w ith the 
optical fluxes, see Fig. \T\ As found by ICzernv et alj <2003fl . the 
medium-energy X-ray flux of NGC 4151 varies on time scales ~ 5- 
10 3 d, whereas the optical variability is also present on longer time 
scales. For hard X-rays, a correlation with the optical and the 1.5- 
12 keV fluxes is less clear. The main peak of the optical emission at 
MJD =! 49000-51000 (1993-1998) is reflected in the 1.5-12 keV 
flux but not in the 20-100 keV one (including the BATSE data). The 
later data from BeppoSAX, INTEGRAL and Suzaku show an overall 
agreement with the optical and softer X-ray fluxes in a sense that 
the minima at MJD -52000 and MJD -54000 and the maximum at 
MJD -52800 appear for all these bands. However, the scarce hard 
X-ray coverage prevents unambiguous conclusions. 

The results shown in Fig. [TJa— b) show that the hard X-ray 
emission of NGC 4151 has varied in well defined limits over the 
last 39 years. Almost all fluxes are within the extreme 20-100 keV 
fluxes from ISGRI (with the maximum in 2003 May and the mini- 
mum in 2007 May). There are four flux measurements from older 
(1975-77) observations above this maximum. However, they are 
relatively doubtful, corresponding to marginal (as 2-3cr) detections 
at > 40 keV. Thus, ISGRI 20-100 keV flux range appears very close 
to the corresponding overall actual range. This is also supported by 
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Figure 2. Medium and hard X-ray light curves from INTEGRALflSGRl 
compared to those from RXTE and Swift (shown with black circles and error 
bars). The 18-50 keV ISGRI flux is shown by open circles on all panels. The 
colours identify the three main states, bright (blue), medium (green), dim 
(red) and, additionally, data excluded from the spectral analysis, Revs. 0076 
and 0809-0811 (magenta), and the rest (cyan), (a) The RXTE/ ASM count 
rate in 30-d bins (right axis), (b) The RXTE/PCU2 count rate (right axis). 
The vertical lines delineate the two periods for which the dim state PCA 
spectrum was obtained, (c) Comparison with the Swift/BAT 15-50 keV flux 
in 7-d bins. The vertical lines show the times of the X4 (dim), X5 (medium) 
XMM-Newton observations and that of Suzaku (dim). The horizontal dashed 
lines show the average flux levels of ISGRI bright, medium and dim states. 



the BATSE data, Fig. |TJb), varying within this range. We also no- 
tice that NGC 4151 was rarely seen at low hard X-ray fluxes before 
the INTEGRAL launch. 

Fig.|2]compares ISGRI 18-50 keV fluxes with the contempo- 
raneous data in the medium and hard X-ray bands from RXTE and 
Swift. The ISGRI fluxes are well correlated with those from both 
RXTE and the Sw i/j/BAIq In particular, we see a good correlation 
with the 1.5-12 keV rate in spite of the variable absorption affect- 
ing a lower part of this band. Fig. [2] also shows that a large fraction 
of the INTEGRAL observations happened during periods with very 
low hard X-ray flux. Fig.[2]also identifies the INTEGRAL observa- 
tions used for the spectral sets (TableQJ. 

The BAT light curve, Fig.[2]c), shows that the hard X-ray flux 
varies by a factor of a few on a week-time scale. When the source 
is sufficiently bright, it is possible for ISGRI to monitor hour time- 
scale variability. This was the case for the bright state of 2003 May, 
and the source showed such variability, see Fig. [3] The shown IS- 
GRI and OMC light curves were extracted with a time bin equal 
to a single pointing duration (10 m-2.5 h, typically - 1 h). When 
observed in dimmer states, ISGRI flux remains constant within the 
measurement errors during a given continuous observation. The op- 
tical flux appears constant on hour-day time scales. Still, it varies 
on longer time scales, and it increases by -30 per cent when the 
source is bright in X-rays (2003 May, 2009 June). 

We have extracted the bright (X1-X3) XMM-Newton/EPlC pn 
light curves in several energy bands in 100-s time bins, and com- 
pared with the OMC 100-s light curves. No variability is seen in 
either light curves on time scales < 1 h. The XMM-Newton light 

2 http://swift.gsfc.nasa.gov/docs/swift/results/transients 
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Figure 3. Illustration of short time-scale variability. The blue bars (left axis) 
and the red dots show the 18-40 keV ISGRI flux and the optical OMC flux, 
respectively, during the dedicated INTEGRAL observations. The vertical 
lines and the green arrows in the top panel show the times of the bright- 
state X1-X3 XMM-Newton and 11 RXTE/PCA observations, respectively. 
The data used for the bright INTEGRAL spectrum are those in the top panel 
simultaneous with the OMC observations except those around MJD 52788 
(Rev. 0076), which have markedly lower fluxes. 



curves vary slowly in good agreement with the trends observed for 
ISGRI (and also JEM-X) fluxes. 

The INTEGRAL data are suitable for correlation studies be- 
cause they provide simultaneous data in the optical, medium and 
hard X-ray bands. Here, we use all the observations with the off- 
axis angle < 3° for which OMC and IEM-X data are available. 
They include also part of the observations not used in our spec- 
tral analysis, namely those identified with magenta circles in Fig. 
E] Fig. [4] shows the correlations between the 18^-0 keV flux and 
those fluxes in the 3-10 keV and V bands. Both show strong corre- 
lations, but the correlation patterns are different. The optical fluxes 
fall into two separate regions, within each they are correlated with 
the hard X-rays. Significant rank-order probabilities (> 0.999) of 
the ISGRI/OMC correlation are found separately for each of the re- 
gions, but also and for the entire data set. The observed bimodality 
may result from strong optical long-term variability between the X- 
ray observations. On the other hand, the 3-10 keV and 18^40 keV 
fluxes exhibit a single, extremely strong, linear correlation. Statisti- 
cal uncertainty blurs somewhat the correlation for low fluxes. This 
result indicates that during all INTEGRAL observations the slope of 
the X-ray spectrum of NGC 4151 below 40 keV remains approxi- 
mately constant. 
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Figure 4. (a) The OMC 5500 A flux vs. the ISGRI 1840 keV flux, (b) The 
JEM-X 3-10 keV flux vs. the ISGRI flux. Each cross corresponds to a sin- 
gle INTEGRAL pointing. The colours identify the three main states, bright 
(blue), medium (green), dim (red), and, additionally, Revs. 0076, 0809- 
0811 (magenta). 



4 SPECTRAL ANALYSIS 

4.1 Assumptions and selection of spectral sets 

Spectral fitting is performed with xspec 11.3 l lArnaudlll996t) . Er- 
rors are given for 90 per cent confidence level for a single pa- 
rameter, Ax 2 = 2.7. The luminosity and distance, D = 13.2 
Mpc, are for Hp = 75 km s" 1 Mpc" 1 . We use the elemen- 
tal abundances of Anders & Ebihara ( 1982) and the photoelectric 
absorp tion cross-sections from B alucinska-Church & McCammorJ 
dl992T) . The Galactic column density is set to N° = 2.1xl0 20 cm" 2 . 
The assumed inclination angle of the reflector and the hot plasma 
(for anisotropic geometries) is i = 45° JDas et al.l20 05). The reflec- 
tor is assumed to be neutral, and relativistic broadening of reflected 
spectra is neglected. Still, in some cases we test for the effects of 
either a variable Fe abundance of the intrinsic absorber and the re- 
flector, a variable ;', or an ionization of the reflector. 

The main assumption used in our selection of spectral sets is 
that the state of the emitting hot plasma is determined by the hard 
X-ray flux. This is justified by our finding that the shape of ISGRI 
spectra for a given revolution is almost constant at a given hard X- 
ray flux. As an additional constraint, we do not split INTEGRAL 
data from a single revolution into different spectral sets. Based on 
this criterion, we select three average ISGRI spectra, as specified in 
Table[T|and Fig. [2] These spectra are accompanied by spectra taken 
with the other instruments. The state is assigned to each of them 
on the basis of the light curves presented in Figs.[2}{3] In particular, 
the high quality of the SwiftfBAT data allows us to determine the 
hard X-ray flux for periods with no INTEGRAL observations, and 
thus to assign the correct state to the existing lower-energy data 
sets (Table |2J. On the other hand, the data from RXTE determine 
by themselves the hard X-ray state. 

The ISGRI bright (B) and dim (D) spectra are of relatively 
good quality, and correspond to the extrema of the measured hard 
X-ray flux. For the flux range in-between, we identify two flux 
ranges. However, we had no corresponding low-energy data sets 
for the upper of them. Thus, we use the lower of them, denoted as 
medium (M). The fitted data sets are as follows. 

The bright state, B. We use the data from 2003 May made 
at low off-axis angles, as shown in the top panel of Fig. [3] The 
spectrum includes a part of Rev. 0074 (without its beginning) and 
Rev. 0075. For these data, we have simultaneous RXTE and XMM- 
Newton observations, see Table [2] The fitted spectral set consists 
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of ISGRI spectrum, the X1+X2 and X3 XMM-Newton spectra, and 
the RXTE PCU2 spectrum summed over 1 1 observations. 

The medium state, M. We sum ISGRI data from Revs. 0634, 
0636, 0678-0679, and use the summed JEM-X 1 spectrum from 
Revs. 0634 and 0636, and the X5 XMM-Newton spectrum. 

The dim state, D. The summed ISGRI spectrum is from Revs. 
0310-0563. We use the RXTE PCU2, the X4 XMM-Newton EPIC 
pn, and the Suzaku XIS0,3 and XIS 1 spectra. 

We have also investigated the effect of including other avail- 
able spectra for each of the data set, whose resulting sets we denote 
'extended'. The extended set B also includes the spectra from the 
INTEGRAL/]EM-X 2 and SPI, and the RXTE/HEXTE clusters 0, 1 . 
The extended set M also includes the SPI spectrum. The extended 
data set D also includes two INTEGRAL JEM-X 1 summed spectra, 
from Revs. 0521-0522 and 0561-0563, the Suzaku PIN spectrum, 
the RXTE HEXTE cluster 1 spectrum, and the corresponding SPI 
spectrum. We discuss the effect of using the extended data sets in 
Section |4. 3 I below. 

In order to limit the complexity of the fitted models, we have 
decided not to use the X-ray spectra at E < 2.5 keV. A two- 
component partial absorption was then sufficient to model absorp- 
tion of the continuum. Using the X-ray spectra in the 0.1-2.5 keV 
band would have required to add a number of spectral lines as well 
as a soft excess component. As we have tested, including that low- 
energy band would not affect our determination of the parameters 
of the hot plasma and of Compton reflection. 

Thus, the EPIC pn, XIS0+XIS3, and XIS1 spectra were used 
in the 2.5-11.3 keV, 2.5-9.5 keV and 2.5-9.0 keV bands, respec- 
tively. The PCA/PCU2 and JEM-X spectra were fitted in the 3-16 
keV and 4—19 keV bands, respectively. The HXD/PIN spectrum 
was used in the 19-60 keV band because the <19-keV data were 
showing a strong excess. For ISGRI, data below 20 keV were al- 
ways excluded, whereas the high-energy limit was in the 180-200 
keV range, depending on our 3cr detection threshold used to se- 
lect good spectral channels. The ISGRI spectra in the staring mode 
from Revs. 0074—75 show an excess of ~10 per cent below 23 keV, 
whose origin remains unknown. It could be related to a specific set- 
ting of the low-energy threshold in the early period of the mission. 
Thus, we used only the data at > 23 keV in this case. The HEXTE 
spectra were used in the 13-160 keV band for the bright state and 
in 23-1 10 keV for the dim state. The SPI spectra were fitted in the 
23-200 keV band. 

While fitting a given spectral set, the parameters of intrin- 
sic continua are the same for all included spectral data, allowing 
only the normalization of the model to vary between them. This 
reflects our assumption that the hard X-ray flux determines the in- 
trinsic state of the source, even for non-simultaneous data. On the 
other hand, absorpti on and the Fe K emission may vary quickly 
jPuccetti et alJI 20071) . and we cannot expect their parameters to be 
the same for different low-energy spectra associated with the same 
spectral set. Thus, we have allowed them to vary between different 
spectra. Our assumed absorber consists of one fully covering the 
source with the column density of N' H , and the ionization parame- 
ter, £ = 4nFJn, where F x is the 5 eV-20 keV irradiating flux and ;; 
is the density of the reflector. For the low-resolution detectors, PCA 
and JEM-X, we assumed f = 0. For the high-resolution detectors, 
EPIC pn and XIS, we allowed £ > 0, and, in addition, we included 
partial-covering neutral absorption, with the column density of N^ 

covering a fraction of f p of the flux. For the ionized absorption, we 

I 1 J k 

use the model by Done et al. ( 1992). Although it is not applicable 

to highly ionized plasma, it is sufficiently accurate at low/moderate 
ionization, such as that found in our data. Absorption for ISGRI, 



Table 3. Spectral fitting results for the e-folded power-law model includ- 
ing reflection. The xspec model is constant*wabs*absori*zpcfabs(pexrav 
+zgauss), and the normalization, K, is given by the photon flux at 1 keV in 



State £ c [keV] 



* 2 /d.o.f. 



264+ 



Bright 
Medium 
Dim > 1325 



+4S 

26 

> 1025 



1.71+ 
1.81+ 



1.81+ 



0.45 

0.6 

0.92 



+0.12 

+ -8.F 
-Ah 

-0.05 


8.7+0-' 


3179/3414 


4.6+ 14 
2 24+8'Sl 
z - -0.02 


1554/1599 
3338/3384 



HEXTE, SPI and PIN spectra is only marginally important, and we 
thus used only the neutral fully-covering N' H at the values fitted to 
the corresponding PCA or JEM-X spectra. The Fe Ka emission is 
modelled by a Gaussian line, with the centre energy, width, photon 
flux and equivalent width of £p e , o~p t , /p e and wf c , respectively. Fit- 
ting the dim state spectra, we have found that adding an Fe K/3 line 
(at ~l.l keV) was significantly improving the fit for the EPIC pn 
and XIS spectra. 



4.2 The cut-off power-law model 

We first fit an e-folded power-law model including a Compton re- 
flection component. This allows us to compare our results to other 
published ones. We use the pexrav model ( Magd ziarz & Zdziarskil 
1995), and our other assumptions are specified above. The results 
are shown in Table[3] where we do not show the fitted absorber and 
line parameters since they are similar to those for our main model 
of thermal Comptonization, see Section l4~3l below. 

Our results regarding T and R are similar to those of earlier 
works. However, a significant difference is found for the e-folding 
energy, £ c , which we find to be high, and compatible with no cut- 
off for two of our states, see Table[3] A finite £ c is found only for the 
bright state. Even then, we obtain values well above those reported 
befor e, which are 30-45 keV (P01), 80-1 80 keV JDe Rosa et al] 
120071) and -210 keV (Z02). As discussed in lZdziarski et al]J2003h . 
an exponential cut-off is much shallower than that characteristic of 
thermal Comptonization, and not sharp enough to model the spec- 
tral high-energy cut-offs observed in Seyferts, in particular in NGC 
4151. Thus, our high values of E c result from the good quality of 
the data below the cut-off, which dominate the statistics and force 
an approximately straight power law to continue to just below the 
beginning of the cut-off. 



4.3 Comptonization model 

To model thermal Compton ization, we use the compps model 
( IPoutanen & Svenssonlll996l) in xspec It models Compton scat- 
tering in a plasma cloud of given temperature, T c , and Thom- 
son optical depth, r, for a number of geometries and locations 
of the seed photon sources. Since T e and r are strongly intrin- 
sically anticorrelated, we also use the Compton parameter (e.g., 
iRvbicki & Lightm an 1979), y = 4(&7' e /m c c 2 )T (where m e is the 
electron mass) as the second fitting parameter instead of t. The 
advantage of this choice is that T t and y are almost orthogonal; 
y determines closely the power-law slope at low energies whereas 
kT s determines the position of the high-energy cut-off. The seed 
photons are assumed to have a disc blackbody distribution (diskbb 
in xspec iMitsuda et al]|2004l) with the maximum temperature of 
kTbb = 10 eV. The model normalization, K, is that of the seed disc 
blackbody model as defined in xspec, see Section l5~3l 
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Table 4. (a) Spectral fitting results with thermal Comptonization as the main continuum. The xspec model is constant*wabs*absori*zpcfabs*(compps 
+zgauss). See Section 1431 for the definitions of the parameters. Lx r and L' x are the bolometric model luminosity (in units of 10 43 erg s ) with and 
without the reflection component, respectively, and only the values for the geometry are given as the values for other geometries are similar. 



State 



kT s [keV] 



R 



K[\0 S 



L%y 



"X7 



Geometry 



r7d.o.f. 



Bright 



54+7 

62+' 

62+' 

73+' 4 
73+13 

73+" 

73+> 5 
61!? 
57+1 3 
57+1° 



1.10+ 



1+0.02 

-0.02 

60+ 001 

61+°-° 4 

-QQ1 

1 06+™ 

1 06+° ° 3 
-Q-PJ 

1 71+0.02 

-0.02 
1 91+0.03 

1,x -Q.Q2 

0.80^ 

0.86+™! 

0.86+!!!!! 



2 6+ 03 

z -°-0.3 

1.3+01 
! 3+8:] 

1 q+8:5 

1,y -Q.2 

1.9+0-3 

9 1+8:4 

Z -0.6 

2 l+°- 3 
z -'-0.4 

1 7+ ai 

'■'-0.2 

1 9+O.2 

' -Q.4 

1 9+0.2 



40+" "' 

-Q.Q6 
55+™ 

-P.-PJ 
55+0-* 

-Q.Q6 
47+0-05 
u '^'-0.07 

0.48+lf 

50+° ° 8 
u - JU -0.06 
n qo+0.04 

039+O.og 
0.38+H1 



5 2+ 1 - 4 

3.9+ 
12.0+ 



5.18 4.61 



-1.6 

5.4+}, 



11.9 

4.7 4 



-1.4 

1.5 
1.0 

1 2+ 01 
6.4+ 06 
4 9+1:5 



sphere (0) 

slab (1) 

slab (-1) 

cylinder (2) 

cylinder (—2) 

hemisphere (3) 

hemisphere (—3) 

sphere (4) 

sphere (-4) 

sphere (—5) 



3165/3414 
3166/3414 
3166/3414 
3168/3414 
3168/3414 
3168/3414 
3168/3414 
3172/3414 
3169/3414 
3166/3414 



Medium 



128+ 



1.02 



-0.13 



1.0 



+0.5 
-0.3 



<0.41 



3.9 



-1.0 



2.26 2.16 



sphere (0) 



1554/1597 



Dim 



190+J 3 
200+] 
188+6 

210+' 2 

215+ 12 

997+10 

-M 

220+ 1 3 

186+^' 
191 + 13 

"'-17 

196+J 3 



0.98+™! 
0.49+™! 
0.48^; 

1 05+™' 

•"-Q.pl 

1 01+™' 

i.ui_ 001 
j 21+001 
1 14+8:81 

0.66+™| 
0.86!;! 
0.8l4°! 



66+ 004 

-Q.Q5 

31+0-01 

UJ1 -0.01 

33+ 001 
-Q.Ql 

64 +tl04 
-Q.Q4 

0.60+g 

0.68+™ 4 
66+'' 03 

u " D -0.04 

45+0-02 

" -0.03 
a SR+0-05 

53+ - 04 



75+0- 04 
"■'-'-0.04 
1,01+005 

, 4+8:8? 

79+0:84 

o: 8 6+|| 

0.84+; 

90 +b - ig 
-BQ5 

77 +«.8+ 

"■"-0.Q4 
q 77 +0.04 
"■"-0.Q4 
77 +0.04 



1 =14+0.01 
,J -Q.Ql 

2 04+ 001 

z -"+-0.02 

! 20+001 
9 o<;+8:82 

-fl-fil 

1 7^+0.01 

•' -Q.Ql 

2 89+ 001 

-8-Q2 

1.71!}! 

0.866+™« 4 
1 56+ 001 

,J -Q.Ql 

1.49+0 ° 



1.47 1.19 



sphere (0) 

slab (1) 

slab (-1) 

cylinder (2) 

cylinder (—2) 

hemisphere (3) 

hemisphere (—3) 

sphere (4) 

sphere (-4) 

sphere (—5) 



3337/3384 
3335/3384 
3350/3384 
3334/3384 
3336/3384 
3334/3384 
3336/3384 
3341/3384 
3337/3384 
3337/3384 



Table 4. (b) Continuation of Table [4ja). Fit results with the thermal Comptonization model (for a sphere, the geometry parameter = 0) regarding the 
absorber, line component, and the relative normalization, Cisgri. of a given X-ray spectrum with respect to ISGRI. See Section l4~T1 for the definitions of 



the other parameters. The units of /Vh, f, and 7p e are 10 22 cm 2 , 10 2 erg cm s 



, respectively, and f ' denotes a fixed parameter. 



State 



X-ray spectrum 



Cisgri 



i : 



/V p 



£ F e [keV] o-pe [keV] 



w Fe [eV] 



Bright EPICpn(Xl+X2) 0.92+ 

EPICpn(X3) 1.08+ 

PCAPCU2(B) 1.17+ 



4 n+u+j 

3.0+[| 
7.4+«- 4 



10+ 



o+- 

64+ 4 ] 

Of 



15 
9.5 



,+3.9 
'-3.3 
+4.0 
-1.1 

Of 



0.39; 

0.55+ 



6.400+ 



1+0.011 
-0.011 
6 393+0.013 

U - J "-0.012 

6 06+ 008 



0.07+ 



7 +0.02 

-0.01 

07+ 001 
"•"'-0.01 

36+° 15 

U.JD 013 



3 0+ 03 
6.8+>; 3 



83+ 



76+' 
169+ 32 



Medium 



EPIC pn (X5) 
JEM-X (M) 



+0.08 
0.06 
-74+O.O6 



83+ 



7 6+ u 3 

13.0+37, 



4+11 

^-3 

Of 



22.6+| 4 
Of ' 



0.50 



+0.05 
0.04 



6.388+ 



0.06 



+0.01 



2.7 



+0.5 



149 



+27 
12 



Dim 



EPIC pn (X4) 
XIS0.3 (D) 
XISl (D) 


1 17+0.01 
1 11+ 001 

111 -0.01 
1 io+0.01 

16+8:81 

llu -0.01 


636 +o.i4 4+30 23 5 +o.6 06 9+o.oi 

"■■'"-0.15 -3 ZJ -0.7 ""'-COl 
m+012 n -7+0.9 90 0+0.9 n <;n+o.oi 
8.03 _ 01fi O.7_ 04 28.8 10 0.50_ 0m 

fitted together with XIS0,3 


fi ano+0.004 
°- 3yS -0.005 
c ^Q/i+0.003 
°- 384 -0.003 

6 412+000 5 

°- H1 -0.004 


06+ 001 

"-OO-0.01 

03+ 001 

"■" J -0.01 

04+ 001 
"■"^-0.01 


PCA PCU2 (D) 


13.6+5j| Of Of 


6.4f 


Of 



2.3+ 
2.2+ 
2.3 
2.1 



+8:3 



271 + || 

291+| 3 
280+] 3 
303+ 3 ^ 



Among the geometries included in compps, we consider here 
four spherical cases; an approximate treatment of radiative transfer 
using escape probability (which is denoted in compps by the geom- 
etry parameter = 0), a sphere with central soft photons (geome- 
try 4), homogeneously distributed seed photons (geometry -4) and 
seed photons distributed according to the diffusion-equation eigen- 
function, oc sin(7TT'/T)/(7TT'/7"), where < t' < t (-5). Then, we 
consider a slab with the seed photons either at its bottom (geometry 
1) or distributed homogeneously (geometry -1). Also, we consider 
the hot plasma in the shape of either a cylinder with the height equal 
to its radius, or a hemisphere, with the seed photons being either at 
its bottom (2, 3, respectively) or homogeneously distributed (-2, 
-3, respectively). 



about Xy emission from brightest Seyferts can be achieved with the 
current satellites. Therefore, we extensively test various variants of 
the Comptonization model for our two best-quality states, bright 
and dim. This could, in principle, give us some indications regard- 
ing the actual geometry of the source. Our results are presented in 
Tables |4ja) and (b), regarding the parameters of the Comptonizing 
plasma and the strength of reflection (which are the same for all 
used detectors within a given state), and of the absorber and the Fe 
line (which are specific to a given X-ray detector within each of the 
states), respectively. The results in Table|4tb) and the entries for the 
bolometric model luminosity, Lx y , in Table|4fa) are given only for 
the case of spherical geometry calculated using escape probability 
formalism (geometry 0). 



The spectral data of NGC 4151 studied in this paper are of 
high quality and represent the largest set ever collected for this 
source. This gives us a possibility to test how precise information 



We find the thermal Comptonization model to provide very 
good fits to the data. The values of the reduced,^ 2 are < 1, which, 
however, appears not to be due to an overly complex model. In par- 
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Figure 5. The 90 per cent confidence regions for the Compton y parameter 
(left) and the optical depth, r (right), vs. the temperature, &r c for three ge- 
ometries of the Comptonizing plasma. The solid blue, solid red and green 
dashed (only for y) contours correspond to the bright, dim and medium 
state, respectively. The dashed curves for r vs. kT e correspond to r oc kT^, 
normalized to the bright-state results. 



ticular, the high-energy continua are determined only by three pa- 
rameters, y, R, kT c . For the bright state, Comptonization provides a 
much better fit than the e-folded power law, with A^ 2 of up to - 14. 
However, any trends seen for the geometry of the Comptonizing 
plasma are rather weak. In the case of the bright state, equally good 
models can be obtained with either a slab or a sphere. Also, models 
with seed photons being distributed within the source are of simi- 
lar quality as those with localized seed photons. On the other hand, 
the differences between the models are slightly stronger in the case 
of the dim state. In particular, the slab model with seed photons at 
its bottom is somewhat better than that with the seed photons dis- 
tributed throughout it. This can be a hint that the actual source has 
the seed photons external to the source rather than internal (which 
would be the case for the dominant synchrotron seed photons). We 
also note that the fitted values of kT c are rather insensitive to the 
assumed geometry in both the bright and dim states. Our results 
for the bright state are similar to those of previous Comptonization 
fits, in particular to those of Z02, who obtained kT^ = 73*™ keV, 
y = 0.88!°! 2 , and R = 0.60!°; 24 at L Xy * 3.9 x 10 43 erg s" 1 using a 
broad-band X-ray spectrum from ASCA and OSSE. 

In all geometries, we see in Table |4{a) that the Compton y 
parameter remains approximately constant between the bright and 
dim states, with only a slight decrease in the dim state. For the 
medium state, Table |4ja) gives the results only for the model with 
escape probabilities, for which the best fit y is almost the same as 
that in the dim state. These results are illustrated by the confidence 
contours shown in Fig. [5] 




50 100 

Energy [keV] 

Figure 6. The Comptonization model spectra (solid curves; the geometry 
parameter = 0) of the three states of the NGC 4151, shown together with 
ISGRI spectra. The bright, medium and dim spectra are shown by the blue 
crosses, green circles and red triangles, respectively. The dashed curves 
show the reflection model components for the bright (upper) and dim state. 
The fitted JEM-X, RXTE, XMM-Newton and Suzaku spectra are not shown 
for clarity. 



Apart from the results given in Table [4] a reduction of x 2 
appears when the Fe abundance (assumed equal for the absorber 
and the reflector) is a free paramet er. We found a moderate over- 
abundance, 1.5+0.2 times that of I Anders & Ebihard dl982h for 
the bright state, with A^- 2 = -6, and 1.3 ± 0.1 for the dim state, 
A^- 2 = -2. For the XMM-Newton observatio n in 2000 December , 
an Fe overabundance by 2-3 was claimed bv lSchurch et alj (2003) 
but that fit was done with the XMM-Newton spectrum only and the 
reflection strength (affecting that determination) was not well con- 
strained. 

We also find that both the bright and dim state spectra prefer an 
inclination angle > 45°, with A^- 2 = -3 and -6, respectively, at 75°. 
Since this effect is observed for symmetric (sphere) and asymmetric 
(hemisphere) geometries, we conclude that its main cause is the 
changing shape of the reflection component. On the other hand, 
allowing the reflector to be ionized does not improve the fit, and 
<f < 0.2 erg cm s _I (at an assumed reflector temperature of 10 5 K). 

The lack of the far-UV and very soft X-ray spectra does not 
allow us to fit the maximum temperature of the seed disc blackbody 
photons, kT\,h- We have found that changing it to either 5 or 20 eV 
yields no improvement of the fit. This is indeed expected given that 
the photon energies emitted by the disc are much below the fitted 
energy range, at which the shape of the Comptonization spectrum 
has already achieved a power-law shape independent of kT bb . 

Table|4jb) also shows the normalization of a given X-ray spec- 
trum relative to that from ISGRI. For the bright state, when the IN- 
TEGRAL, XMM-Newton and RXTE observations were almost con- 
temporary, Cisgri is ver y close to 1 and, for both of the XMM- 
Newton spectra, it follows the flux changes seen in the top panel 
of Fig. [3] This shows that ISGRI, EPIC pn and PC A detectors are 
well cross-calibrated. In the case of the two other spectral sets, the 
observations are not simultaneous. Nevertheless, Cisgri varies only 
within the range of 0.74-1. 17, confirming the validity of our selec- 
tion of the spectra for a given flux state. 

We have also tested our results using the extended data sets 
(see Section |4~TT i. However, the quality of the fits and the plasma 
parameter uncertainty are not improved when they are used. This 
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Figure 7. The Comptonization models for the extreme states (solid curves, 
same as in Fig. [6) compared to the average OSSE spectrum (crosses and two 
2<x upper limits). The rightmost upper limit is the COMPTEL upper limit in 
the 0.75-1 MeV band ( Maisack et al. 1995). The histograms give the INTE- 
GRAL/PICsYT 3<x upper limits for 1 Ms (upper) and 10 Ms (lower) effective 
exposure time. The PICsIT 1.34 Ms average 277^4-61 keV flux is shown 
by erro r bars with a diamond. The SPI spectrum obtained bv lBouchet et alj 
12008) is shown by error bars with circles. Note that the 200-600 keV data 
point may be spurious. 



confirms that our selection of the primary sets was valid. Thus, we 
do not present here those results. 

We note that even for such a bright AGN as NGC 4151, we 
have no data at energies > 200 keV, as shown in Fig. [6] Still, 
for the dim state, we have found kT e ^ 200 keV with a statis- 
tical uncertainty of ~ 7 per cent or so. The physical reason for 
this low uncertainty is the observed lack of spectral curvature. 
In the process of thermal Comptonization, a power law slope 
at E < kTe appears due to the merging of many individual up- 
scattering profiles, each being curved. The cut-off at E > kT e 
is due to the electron energies being comparable to the photon 
ones, in which case a photon no longer increases its energy in a 
scattering. Since we see no cut-off up to ~200 keV, models with 
a lower kT e are ruled out. On the other hand, models with a 
higher kT c correspond to an increase of the photon energy in a 
single scattering by a factor so large that the individual, curved, 
scattering profiles become clearly visible and the spectrum at 
E < kT c is no longer a power law. Since the observed spec- 
tral shape is a power law, this case is also excluded. We note, 
however, that our assumed model corresponds to a single scat- 
tering zone. A power law spectrum at E < 200 keV may also be 
produced by a superposition of emission regions with different 
temperatures (smoothing out the resulting spectrum), as, e.g., 
in a hot accretion flow. In this case, the maximum temperature 
of the flow could be higher than 200 keV. In any case, we do 
see a very clear difference in the characteristic Comptonization 
temperature between the bright and dim state. 

We comment here on the very high value of kT s = 315 ± 15 
keV obtained by P01 for the 1999 January BeppoSAX spectrum of 
NGC 4151. The fit of that model appears rather poor when com- 
pared with the curvature of the BeppoSAXfPDS data above 50 keV, 
see fig. 2(b) in POL We also find tha t the PDS spectrum is very 
similar to the average OSSE spectrum {Johnson et alj 19971) , shown 
in Fig. [7j (with only the normalization of the PDS spectrum being 
about 5 per cent lower). This indicates that the source was in a mod- 



erately bright state in 1999 January, for which state the temperature 
of 315 keV appears highly unlikely. We cannot explain this dis- 
crepancy; it may be that kT c = 315 ± 15 keV represented a spurious 
local minimum. 

As we see in Fig. [7] our bright-state model sp ectrum agrees 
well in shape with the average OSSE spectrum fljohnson et alj 
1 19971) . The latter also agrees with the first three po ints of the av- 
erage NGC 4151 spectrum from INTEGRAL/SPl JBouchet et ail 
l2008h . However, the SPI 200-600 keV flux is implausibly high, 
which may reflect a problem wi t h the b ackground estimate. Indeed, 
the SPI analysis of IPetrv et alj ([2009) yields only upper limits at 
E > 150 keV. 

We see in Fig.|7]that the Comptonization models of the bright 
and dim states cross at = 300 keV. Unfortunately, it is impossible 
to verify it with our data. In particular, we cannot test if there is any 
non-thermal component at high energy. Even the average spectrum 
from OSSE, the most sensitive detector up to date at E > 150 keV, 
gives only a weak detection of NGC 4151 at £ > 250 keV with 
an > 4 Ms exposure time. From the shown upper limits (including 
systematic errors) of the INTEGRALfPlCsYT, the current detector 
most sensitive at E > 250 keV (Lubinski 2009), we see it would 
need > 10 Ms exposure time to detect NGC 4151 in the 300^100 
keV range. Merging all the suitable INTEGRAL data (with the ef- 
fective exposure time of 1.34 Ms taken without the bright state, for 
which the observation was in the staring mode, not suitable for the 
PICsIT), we extracted the PICsIT flux in the 277-461 keV band. 
Although it is not a detection (and the shown error bar does not 
include systematic errors), it provides a hint that the average NGC 
4151 emission in this energy range is weak. 



5 DISCUSSION 

5.1 Reflection and Fe line 

5.1.1 Comparison with previous work 

A reflection compone nt in the X -ray spectra of NGC 4151 was 
first found by Zd ziarski et alj dl996af) . using contemporaneous data 
from CGRO and Ginga, yielding R ^ 0.43 (using a Comptoniza- 
tion model and assuming i = 65°). Then, Z02 found R = 0.6 (at 
i = 17°) using CGRO and ASCA data, and P01 found R a 0.2- 
1.8 using BeppoSAX data and a reflection model averaged over 
the viewing angle. A previous analysis of INTEGRAL bright-state 
data from 2003 May y ielded a larger reflection, R = 0.72 ± 0.14 
JBeckmann et alj |2005). than that determined by us (their value of 
kT c = 94+ 4 keV is also different from that found by us for the 
same data). By fitting the ISGRI and JEM-X bright-state spectra 
alone, we find results very similar to those with the XMM-Newton 
and RXTE spectra included (Tabled- Therefore the differences ap- 
pear to be a consequence of the substantial change o f the INTE- 
GRAL calibration since 2005. Furthermore, Beckmann et al. (2005) 
used the SPI spectra obtained in the staring observation mode, 
which can be affecte d by a large background uncertainty. Then, 
|Pe Rosa et alj J2007I) found R ^ 0.4 for several BeppoSAX spectra 
and R =* for the remaining ones, but usin g an e-fol ded power- 
law model rather than Comptonization. Then. lSchurch & Warwick! 
(2002) found R ^ 1.9 using XMM-Newton data assuming i = 65°, 
but this result was obtained using the 2.5-12 keV band only and 

assuming a fixed T = 1.65. ^_^_^ 

Regarding Fe Ka line, Izdziarski et alj ll996ah found a nar- 
row line with / Fc * (0.6-1.0) x 10~ 4 
and Z02 found a narrow line at /p c - 2.5 x 10~ 4 cirT 
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companied by a weaker broad relativistic line using ASCA data. 
The broad component has not been observe d later, whereas the nar- 
row li ne flux measure d by BeppoSAX (De Rosa et al. 2007), Chan- 
dra dOgle et alfeOOd) and XMM-Newton dSchurch et al.l2003h was 
found to be in the range (1.3^4.2) x 10~ 4 cirT 2 s . As seen in Ta- 
bles Efa-b), our results on R and / Fi . approximately agree with the 
previous measurements. 

5.1.2 Close and distant reflector 

We see in Table |4{a) that the relative strength of reflection, R, is 
close to being twice higher in the dim state than in the bright one. 
(In the medium state, the amount of data is very limited, and the 
apparent lack of detectable reflection may be not typical to that 
flux level.) The increase of R with decreasing flux could be a real 
change of the solid angle subtended by the reflector between the 
states. On the other hand, it could be that it is due to the contribution 
of a distant reflector, in particular a molecular torus, as postulated 
in the AGN unified model. Reflection from distant m edia is indeed 
commonly seen in Seyferts 2, e.g., iRevnolds et al.l dl994l) . If we 
attribute the increase of R in the dim state due to that component, 
we can calculate the fractional reflection strength from the close 
reflector, presumably an accretion disc, and the flux reflected from 
the torus, which, due to its large size, is averaged over a long time 
scale, and thus assumed to be the same in either state. 

The hypothesis of the stable geometry of the source is sup- 
ported by the approximate constancy of the Compton y parameter, 
see Table [4ja) and Fig. [5] and of the photon index, T, see Table 
[3] Those parameters are closely related to t he amplification factor 
of Comptonization (e,g.. lBeloborodovl l999). whose constancy can 
be most readily achieved if the system geomet ry is constant (e.g., 
lHaardt & Maraschill99lllZdziarski et aljl999l) . Thus, we consider 
an approximately constant relative reflection strength from the disc 
very likely. Observationally, we do see a strong R-T correlation in 
Seyferts and X-ray binarie s, which can be int erpreted by a geo- 
metrical feedback model dZdziarski et al.lll999l) , and which would 
require a constant disc R for a given T. 

The unabsorbed flux in a given state, F$ , where S = either D 
(dim) or B (bright), is the sum of the (unabsorbed) incident (i) and 
reflected (r) fluxes, 



Fls + F r „ 



(1) 



where Fs and F ti $ can be calculated from the results of our spectral 
fitting in Section |4~3l We assume that the observed (unabsorbed) 
reflected flux in each state is the sum of the close disc (d) state- 
dependent component and the constant distant torus (t) component, 



F r ,s - F d , S + F t . 



(2) 



Then, assuming the constancy of the solid angle subtended by the 
disc reflector, or, equivalently, its reflection strength, R d , the disc- 
reflected fluxes are, 



F d ,s - F TtS (R A /R s ), 



(3) 



where R$ is the reflection strength for a given state. We have four 
equations {2}{3) for four unknowns, F dB , F dD , R d and F t , which 
can be readily solved. We note that they can be formulated without 
involving Fys , which implies that the solution does not depend on 
the choice of the energy interval in which the fluxes are measured 
as long as it includes the entire reflected spectrum. The solutions 
for Rj and F, are, 

(F,b - F lD )R B R D (R D —Rn)F T $F T p 



R 



F T fiRo - F r pR$ 



F,,bRd — F T pR% 



(4) 



Then, the fractional reflection strength of the torus can be approxi- 
mately estimated as, 



R t - («>(F t /(F r », 



(5) 



where (R) and (F r ) are the average observed reflection strength and 
the average reflected flux, respectively. 

Using the values obtained for thermal Comptonization using 
escape probability formalism (geometry parameter = 0), we ob- 
tain R d = 0.27 ± 0.07, F, = (8.2±0.6)xl0- u erg cnr 2 s _1 , and 
R t = 0.24±0.04 using either arithmetic or geometric averages in 
equation ((5). Similar values are obtained for other Comptonization 
geometries considered in Section POl We can check the robustness 
of our estimates by allowing the reflection strength to depend on the 
accretion rate, or flux. If we assume that R d in the dim state is a half 
of that in the bright state [which requires an appropriate change of 
equation (|3)L we obtain relatively similar values of the disc reflec- 
tion in the bright state of F d , B = 0.23±0.12, F, = (l.l±0.1)xl0-'° 
erg cm" 2 s" 1 , and R t = 0.32±0.04. 

Similar reasoning can be applied to the Fe Ka line. The ob- 
served line photon flux in the state S, Is, is assumed to be the sum 
of contributions from the disc, 7 d> s and the torus, /,. The latter is 
assumed constant, and the former proportional to the ionizing flux 
incident on the disc. As a simplification, we assume the ionizing 
flux (> 7. 1 keV for neutral Fe) to be proportional to the differen- 
tial continuum photon flux, F(E)/E, at the line centroid energy, 
which we denote N$, and which is equal to /s/w s , where ws is 
the observed equivalent width in a given state. This then implies a 
constant disc line equivalent width, w d . Above, for the sake of sim- 
plicity, we have dropped the indices 'Fe', used in Section |4~3l and 
Table|4jb). The equations are, 



4s + 4 



h 



■ w d N & , 



which can be solved for 

h - Id r IdNb 



»'j 



N B -Nu 



/, 



hNn 



Nn-Nv 



(6) 



(7) 



The equivalent width of the torus line flux with respect to the aver- 
age incident photon flux, {N), can be estimated as, 



h/(N). 



(8) 



For the numerical values, we use here averages of the fit re- 
sults for both EPIC pn spectra in the bright state, and for the EPIC 
pn and XIS spectra in the dim state, see Table |4{b). We find 1% = 
(2.9±0.2) xl0~ 4 cm" 2 s _1 , 1 D = (2.3±0.1)xl0~ 4 cnr 2 s _1 , N B = 
(3.6±0.1)xl0~ 3 keV- 1 cnr 2 s" 1 and N D = (8.1±0.1)xl0- 4 keV" 1 
cnr 2 s" 1 . This yields w d = 23±11 eV, /, = (2.1±0.9)xl0~ 4 cnr 2 
s~', and w t = 93±41 eV, 120±53 eV using the arithmetic or geo- 
metric average, respectively. 

The equivalent width of the Fe line at i = 45° by an 
isotropically illuminated cold disc for T = 1,75 is - 140 eV 
(George & Fabian 1991). The value of R d a 0.27 found based on 
the observed reflection strength then implies w d ~ 38 eV, some- 
what higher than the estimate based on the line fluxes, but still in 
an approximate agreement taking into account measurement errors 
and a number of assumptions we have made. On the other hand, 
R t 22 0.24 found above would explain only w, ~ 30 eV. However, 
we note that we have neglected the local absorber, which obviously 
also gives rise to an Fe Ka line component. Its characteristic A'h 
of ~ 10 23 cnr 2 (Table Efc) can readily explain the excess equi va- 
lent width of ~ 70 eV (e.g.. lMakishim.Jl986l : lA"waki et aljl99lh of 
the constant line component with respect to that expected from the 
torus. 
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We note a number of cave ats for our results. The st andard 
accretion disc is flared, not flat dShakura & Sunvaevll 19731) , This, 
however, would have a relatively minor effect, changing some- 
how the distribution of the inclination angles. Our results indi- 
cate disc reflection with Ra substantially less than unity, which 
implies the Xy source is not entirely above the disc. A likely 
geo metry explaining it is a hot inner flow surrounded by a disc 
fe.g. lAbramowicz et alJ[l99l ; lNaravan & Yill 1 9951 ; I Yuarjl200 lh . in 
which case the incident radiation will have much larger incident an- 
gles (measured with respe ct to the axis of symmetry) th an those as- 
sumed in the used model (Magdziarz & Zdziarski 1995). The disc 
may be warped (e.g.. IWiiers & Pringle 1999), which again would 
change the distribution of the incident angles. Furthermore, Comp- 
ton scattering in the ho t plasma above the disc reduces the observed 
reflection strength (Pet rucci et alj|2001bh . 

Then, we also used the slab geometry for the torus reflection. 
Thus, the obtained solid angle does not correspond to the actual 
angle subtended by the to rus from the X-ray source. For example, 
IMurphv & Yaqoobl (2009J) considered a torus with a circular cross 
section subtending (as seen from the centre) a 2n solid angle. They 
found that the reflection component in this case is several times 
weaker than that corresponding to a slab subtending the same angle. 
One obvious effect here is that, unlike the case of a slab, an observer 
sees only a fraction of the reflecting surface, without parts obscured 
by the torus itself. If we take this into account, our value of R t ~ 0.2 
appears consistent with a torus subtending a. ~ 2n solid angle. We 
note, however, that even if the torus solid angle is formally 27t or 
so, a large part of it will be shielded from the X-ray source by 
the accretion disc and the black hole itself, so the actual irradi- 
ated solid angle may be substantially lower. This effect was not 
taken into account in the geometrical model of IMurphv & Yaqoobl 
(2009). Furthermore, the cross section of the torus may be sub- 
stantially different from circular, which may increase the observed 
reflection, as well as it may be cl umpy JKrolik & Begelmanll 19881 : 
iNenkova. Ivezic. & Elitzurl l2002). which will decrease it. Still, our 
value of w t ~ 100 eV agrees with that from a torus with the col- 
umn density of N a ~ 10 24 cirT 2 for T = 1.75 ( IMurphv & Yaqoobl 
2009), which provides an explanation for that equivalent width al- 
ternative, or additional, to that as being due to the line emission of 
the absorber (discussed above). 

We have assumed that the torus reflection co mponent is con- 
stant o n the time scale of years. On the other hand. lMinezaki et al.l 
(2004) found that the dusty torus inner boundary is at =; 0.04 pc. 
Thus, a fraction of the torus reflection may vary on the correspond- 
ing time scale of ~ 50 d. Furthermore, 0.04 pc corresponds to 
= 2 x 10 4 R g (where R g = GM/c 2 is the gravitational radius), where 
the accretion disc may still be present and join onto the torus. On 
the other hand, Radomski et all ((2003) have constrained the torus 
outer boundary to < 35 pc, so the bulk of the torus reflection may 
still be constant over a time scale of years. 



5.2 Absorber properties 

We compare our results with those based on BeppoSAX, which 
also provided broad-band spectra, allowing to simultaneously de- 
ter mine absorber proper ties a nd the continuum. Th ey were studied 
by IPuccetti et ail J2007I) and |Pe Rosa et all J2007I) , who used the 
same absorber model as in this work. We find their results to be 
compatible with ours. The fully covering absorber has Nh = (0.9- 
9.4) x 10 22 cm" 2 for BeppoSAX and =* (3.0-8.1) x 10 22 cirT 2 in our 
case, see Table|4jb). For the partially covering absorber, A^ ^ (3.5- 
30.3) x 10 22 cm- 2 (BeppoSAX) and = (9.5-30.2) x 10 22 cnr 2 (this 



work). The covering fractions are also similar, f v m 0.34-0.71 
(BeppoSAX) and 0.36-0.71 (this work). The only exception is the 
long-exposure BeppoSAX ob servation of NGC 41 51 in December 
2001 showing a very low A^h dDe Rosa et alj 2007). 

We find an anticorrelation between the Nh of both absorber 
components and the hard X-ray flux. We define the total column, 
Nh = N' H + f'N^. For the bright state, we have A'h/IO 22 cnr 2 
= 11.1±3.5 (XMM-Newton X1+X2) and 8.2±4.9 (X3), for the 
medium state, 18.0±2.8 (XMM-Newton X5), and for the dim state, 
23.1+0.8 (XMM-Newton X4) and 23.2±0.9 (Suzaku). Also, the 
covering fraction is anti-correlated with the hard X-ray flux, see 
Tablej4jb). A similar trend is seen for the results of IPuccetti et al.l 
(2007), who used the 6-10 keV for the X-ray flux. The variable part 
of the absorber needs to be relatively close to the X-ray source to 
be able to follow its flux on the time scale of days. The absorber in 
NGC 4151 was identi fied with massive out flows from the accret ion 
disc (P iro et al.ll2005l) . a broad-li ne region l Puccetti et alj|2007h or 
the surface or wind of the torus { [Schurch & Warwickl l2002). Only 
the disc wind is close enough to the X-ray source to follow the 
change of the nuclear emission and to produce any correlation. The 
physical explanation for the correlation remains unclear; possibly 
the wind rate in NGC 4151 decreases with increasing accretion rate. 



5.3 The nature of the X-ray source 

The geometry and parameters of the source are constrained by a 
number of our findings. We find (i) that both the Compton param- 
eter and the X-ray spectral index are approximately the same in 
both bright and and dim states (Sections 14.2444.3b . This implies an 
approximately constant amplification ratio of the Comptonization 
process, i.e., the ratio of the powe r emitted by the pl asma to that 
supplied to it by seed photons (e.g jBeloborodovl 19991) . The ampli- 
fication factors obtained with the Comptonization model are indeed 
similar, A = 17±2, 15±2 for the bright and dim state, respectively 
(assuming kT^ = 10 eV). A similar value of A as 13 was obtained 
for a Comptoniza tion model fitted to the 1991 data from ROSAT, 
Ginga and OSSE dZdziarski et"ai1ll996j) . 

If the seed photons are supplied by an accretion disc surround- 
ing the hot plasma, this implies an approximately constant disc in- 
ner radius. Based on this, we have inferred (ii) a relatively weak 
reflection from the disc, R& a* 0.3 (Section l5.1.2t . This qualitatively 
agrees with A » 1 and both findings rule out the static disc corona 
geometry; the (outer) disc subtends a small solid angle as seen from 
the X-ray source, and the X-ray source subtends a small sol id an- 
gle as seen from inner parts of the disc dZdziarski et alj|l999l) . An 
implication of the latter is that the modelled disc blackbody emis- 
sion (which provides seed photons for Comptonization) should be 
much weaker than that observed (which corresponds to the entire 
disc emission). This seems to be indeed the case as shown in Fig. [8] 
where we see that the bright-state UV flux inferred from the Comp- 
tonization model is about an order of magnitude below the shown 
maximum observed far UV (1350 A; 9.2 eV) flux. On the other 
hand, the model dim-state far UV flux is close to the historical 
minimum observed. However, that minimum repr esented a single 
isolated dip in the light curve dKraemer et al . 2006), and the actual 
far UV flux corresponding to the dim state is likely to be signifi- 
cantly higher. In choosing the shown range of the far UV flux we 
used its strong correlation with the op tical flux, which can be seen 
by comparing Fig.QJd) with fig. 1 in lKraemer et all d2006l) . 

The IR, optical and UV fluxes shown in Fig. [8] appear rela- 
tively weakly affected by the host galaxy emission. In particular, 
the dominance of the AGN in the U band is shown by its strong 
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Figure 8. The broad-band spectrum of NGC 4151. The radio, IR, optical 
and UV data (black dots) were taken from the NED database. Additional 
data are sh own in circles in radio (gree n, Ulvestad et al. 2005), mid and near 
IR (brown. iRadomski et al]2003l ; lRuiz et alJ2003l optical (magenta, OMC, 
this work) and UV (purple, Swj'/f/UVOT, this wor k; and cyan, HST/ST1S, the 
flux extrema from fig. 3 of lKraemer et alj|2006l) . The XMM-Newton EPIC 
pn and ISGRI spectra of the bright and dim states are shown by the blue and 
red dots, respectively, the corresponding unabsorbed Comptonization mod- 
els below 10 keV are shown by dashes, and the disc blackbody (assuming 
geometry parameter = and kT bb = 5 eV) incident on the hot plasma are 
shown by the dotted curves. The upper limits in t he 0.75-1-3 MeV energy 
ranges are from COMPTEL JMaisack et all 19951) . 



variability dCzernv et alj|2003l) . In the optical range, there can be 
some non-negligible fraction of emission from the broad and nar- 
row line regions, but in the UV we expect that the disc dominates. 
Given that the measured UV emission is, furthermore, absorbed by 
the host galaxy, the conclusion above that only a small fraction of 
the disc emission undergoes Comptonization appears secure. 

Another result is (iii) the connection of the normalization of 
the fitted disc blackbody seed emission to the disc inner radius, R m . 
Based on the definition from xspec, we can express it as 

10 8 rcosi(fi in /10 12 cm) 2 



K 



f 4 

■'col 



(O/10Mpc) 2 



(9) 



wher e / co i ~ 1 .7 is the c olour correction to the blackbody tempera- 
ture dShimura & T akahara 1995), and r is the ratio of the blackbody 
emission incident on the plasma to that emitted by the disc. It im- 
plies (at D = 13.2 Mpc and i = 45°), 



Rh 



\10 9 / lO.l/ \ 1.7 



(10) 



At M = 4.6 x 10 7 M o JBentz et alj|2006l) . R g 
Thus, the values of A* in Table[4{a) for the bright state (at T bb = 10 
eV) imply R m ~ 6R, g , the last stable orbit for a non-rotating black 
hole. This is in conflict with the findings above, according to which 
the inner part of the flow is occupied by the hot plasma and not 
by the disc, unless the black hole is rotating fast. Furthermore, the 
diskbb model used for the seed photons becomes invalid when R in is 
close to th e last stable orbit with the standard zero-stress boundary 
condition aGierliriski et alJI 19991) . We note, however, the disc flux 
is oc if? 75,, and approximately R iu oc r^ 2 for a given Xy spectrum. 
In particular, K = 7.2 x 10 9 for the bright state at kT bb = 5 eV 
(consistent with the UV data) and the geometry parameter = 0, 
confirming the above scaling. This yields R m ^ 18_ft g , allowing 
the existence of a hot inner flow. Given that agreement, we used 



kT bb = 5 eV in the models shown in Fig. [8] We note that if R^ ^ 
constant, kT bb has to be somewhat lower in the dim state. 

For comparison, the radius, R CLI , at which the integrated grav- 
itational energy dissipation at R > i? cq equals that at R < R c{] 
is R eq ^ l4.7R g for a n on-rotating black hole in the pseudo- 
Newtonian approximation dPaczvnski & Wi ita 1980) and with the 
zero-stress boundary condition at the last stable orbit, and less for 
a rotating black hole. Thus, R^ n ~ R cq , as found above, is consistent 
with the observed energy distribution of NGC 4151, showing the 
optical/UV emission (presumably due to the disc) at a similar level 
as that at ~ 100 keV (due to the hot flow). Note that R[ n moderately 
larger than R el] is required if some fraction of the energy dissipated 
below R m is advected to the black hole instead of being radiated 
away. However, a disc truncated at R m » R^ would be in strong 
conflict with the broad-band spectrum, as it implies the disc UV 
emission much weaker than that of X-rays. A prediction of this re- 
sult is a presence of an Fe Ka line only moderately relativistically 
broadened (as found earlier by Z02). 

The reflection of R A = 0.3 may be easily achieved if the central 
hot region (with a large scale height) is surrounded by a disc. We 
have calculated the reflection strength in some simple geometrical 
models, with a sphere and cylinder surrounded by a flat disc. For the 
sphere and cylinder with the height equal to 1/2 of its radius we find 
Ra - 0.29 and 0.26, respectively, in agreement with our estimate. 
At this Ri, the emission of the parts of the disc close to the hot flow 
is most likely dominated by reprocessi ng of X-rays, in agree ment 
with the constancy of the v parameter dZdziarski et alJI 19991) . As 
a caveat, we point out that the geometrical R4 estimates are highly 
simplified, not taking into account the actual hot flow geometry and 
the radial distribution of the Comptonized emission. We also note 
that the constancy of the y parameter can also be explained in a 
model in which the inner disc radius increases in the dim state, but 
at the same time the energy dissipation in the hot flow becomes 
relatively stronger in its outer region, e.g., due to advection. 

On the other hand, the effective solid angle subtended by the 
hot flow as seen from the disc can be identified with our param- 
eter r ~ 0. 1 above. All the above considerations very strongly 
argue for the Comptonizing plasma in NGC 4151 being photon 
starved. This is not compa tible with a static corona above the disc 
dHaardt & Maraschilll99ll) . and compatible with a hot inner flow 
surrounded by an outer cold disc. However, it is also possible that 
the corona does cover the disc but it is in a state of a mildly relativis- 
tic ou tflow (Beloborodov 1999; Malz ac, Beloborodov & Poutanenl 
1 200 II) . The relativistic motion then both reduces the flux of the 
up-scattered photons incident on the disc, which explains the rel- 
ative weakness of the disc reflection, and reduces the flux of the 
disc blackbody photons entering the corona in its comoving frame, 
which explains the photon starvation of the hot plasma, i.e., A » 1 . 
A prediction of this model is the presence of a very broad relativis- 
tic Fe line (as the disc now extends all the way to the minimum 
stable orbit). The lack of detection of such a line in high-quality 
XMM-Newton and Suzaku spectra represents an argument against 
this model. 

Our next finding has been that (iv) the plasma temperature 
increased by a factor of =* 3.5 when the bolometric X-ray luminos- 
ity decreased by about the same factor, see Table [4ja). In agree- 
ment with the approximate constancy of y, the optical depth also 
decreased by a similar factor. This can be explained if the accretion 
rate, M, is oc Lx y and the hot flow has a velocity profile independent 
of M, so t oc L Xy - The temperature then follows from the energy 
balance provided the ratio of the power in the hot flow emission to 
that in seed photons cooling remains constant. However, models of 
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accreting hot flows yield a dependence of r on L Xy weaker than the 
si mple proportionality [e.g., eqs. (10-11) in Zdziarski 1998; table 
2 in lYuan & Zdziarskill2004ll . 

The total model luminosity of the Comptonization component 
(including reflection) in the bright and dim state is L Xv — 5.2 x 
10 43 erg s" 1 and 1.5 x 10 43 erg s~ x , respectively (see Tabled). The 
estimated average luminosity from the IR to UV is <* 7 x 10 43 erg 
s _1 , with the IR band (< 10' 4 Hz) contributing ~ 60 per cent of it. If 
we then assume that the IR flux remains constant and the optical- 
UV flux scales as the high-energy emission, the total bolometric 
luminosity can be estimated as 1.3 x 10 44 erg s , 7 x 10 43 erg s 
in the bright and dim state, respectively. Assuming the hydrogen 
mass fraction of X = 0.7 and M = 4.6 x 10 7 Mq, the Eddington 
luminosity, L E = 87tGMm p c/[(l + X)cr T ] * 6.8 x 10 45 erg s~\ 
giving the Eddington ratios of as 0.01-0.02. Note the factor of 2 
variability of the bolometric luminosity compared to the factor of 
3.5 variability of the Comptonization component. This difference 
appears to be explained by the presence of the constant, averaged 
over very long time scales, IR component. 

Relevant time scales include the light crossing time across 
the circumference of the last stable orbit, which for a non-rotating 
black hole is ~ 2700 s (at the default mass). The corresponding 
light crossing time at R m ~ l5R g is ~ 7000 s. Thus, the hour-time 
scale variability observed by ISGRI in the bright state (Fig.[3]( ap- 
pears to correspond to fluctuations at those radii. On the other hand, 
the viscous time of an accretion flow is, 



GM(R\ 312 ,,/tfy 



(11) 



=* 1.2xl0"s 



3/2 



M 



a y 1 (H/R 

\5R.J \ -1.6 x i0 7 M o /VaT/ \~TJ3 



where H is the flow scale height, a is the viscosity parameter, and 
H/R ~ 1/3 for a hot flow. At our default mass, 15/?„ and a = 0.1, 
fvisc - 13 d. For comparison, the Swift/B AT data for NGC 4151 
need to be rebinned to several days in order to get a sufficient 
signal-to-noise ratio for low flux states. Thus, this detector can sam- 
ple the overall behaviour of the hot flow but not variability close to 
the last stable orbit. 

The spectrum of NGC 4151 resembles spectra of black-hole 
binaries in the hard state, also commonly fitted by thermal Comp- 
tonization. Howe ver, a faint non-therma l tail has occasionally bee n 
seen there, e.g.. lMcConnell et al.l d2002l) . IWardziriski et alj J2002h . 
The Comptonization model used by us also allows to include a cor- 
responding non-thermal tail to the electron distribution. However, 
the data for NGC 4151 do not allow us to test for the presence of 
any non-thermal spectral component, see Section POl 



5.4 Comparison with other objects 

We have presented the average parameters of Seyfert spectra in 
Section [T] for both e-folded power law and Comptonization mod- 
els. Our results for NGC 4151 are within the range observed for 
other Seyferts, as noted by Z02. Z02 also provide a critical analysis 
of earlier finding of very hard X-ray spectra in NGC 4151, which 
they find to result from improper modelling of its complex absorp- 
tion. NGC 4151 in the bright state is only slightly harder than the 
average Seyfert 1, and it shows less reflection, in acc ordance with 
the reflection-index correlation dZdziarski et al.l 19991) . On the other 
hand, the dim state appears to have a significant contribution to re- 
flection from a distant reflector, possibly a torus (Section|5TL2]l. 
The similarity between spectra of Seyferts and black-hole 



binari es in the hard spectral state was pointed out by, e.g., 
Izdziarski et alJ dl996b). though the latter appear to have on average 
somewhat harder spectra. This is likely to be due to the disc seed 
photon energies being much higher than those in Seyferts, which 
reduces the spec tral index, T, for a given Comptonization amplifi- 
cation factor, A (Beloborodov 1999). 

A striking similarity between the OSSE spectra of the black- 
hole b inary GX 339-4 in the hard state and NGC 4151 was pointed 
out bv lZdziarski et alj dl998h . lWardzinski et ai] d2002h have shown 
that when L Xy of GX 339^4 decreased by a factor of ~ 2, kT t in- 
creased significantly (see their fig. la), from 46^ keV to 76^ keV, 
with only a slight softening of the X-ray slope. This effect is very 
similar to that seen in NGC 4 1 5 1 , see TableHa). In both NGC 4151 
and GX 339-4, this effect may be due to a decrease of M in the hot 
flow causing the associated reduction of r. 



6 CONCLUSIONS 

We have presented a comprehensive spectral analysis of all INTE- 
GRAL data obtained so far for NGC 4151, together with all con- 
temporaneous data from RXTE, XMM-Newton, Swift and Suzaku 
(Section|2]l. Our main findings are summarized below. 

We have found that the 20-100 keV emission measured by IN- 
TEGRAL has had almost the same range of fluxes as that measured 
by other satellites during past 40 years (Section[3). Thus, our analy- 
sis appears to explore the full range of the variability of this object. 
Also, we have found that this flux appears to uniquely determine 
the intrinsic broad-band spectrum. 

Simultaneous observations by INTEGRAL in the optical range 
and in medium and hard X-rays show a very strong correlation 
within the X-ray band, and a less clear correlation between the op- 
tical and X-ray emission, with a bimodal behaviour of the optical 
flux. The linearity of the former correlation shows that the X-ray 
spectral slope of NGC 4151 remains almost constant. 

Most of the INTEGRAL observations correspond to either a 
bright or dim hard X-ray state (Section [4.U . We have found that a 
thermal Comptonization model provides very good fits to the data 
(Section |4~3l , As the state changes from bright to dim (with L Xy de- 
creasing by a factor of =^ 3.5), the plasma electron temperature in- 
creases from kT t ~ 60 keV to ~ 200 keV and the Thomson optical 
depth decreases from t ~ 2 to ~ 0.6 (for spherical geometry), i.e., 
t oc 1/kTe and t oc L Xy . The former proportionality corresponds 
to an almost constant Compton y parameter, with the X-ray slope 
remaining almost constant, with only a slight softening in the dim 
state, in agreement with the medium/hard X-ray correlation. This 
is suggestive of almost constant source geometry and amplification 
factor. The latter proportionality may occur due to the accretion rate 
varying at an approximately constant hot flow velocity. 

The fitted strength of Compton reflection increased from R = 
0.4 to R - 0.8, which, at the face value, would be in conflict with 
the constant source geometry indicated by the constant y. However, 
in accordance with the AGN unified model, NGC 4151 is likely to 
possess a remote torus, also reflecting the central X-ray emission. 
Assuming the solid angles subtended by both the (close) disc reflec- 
tor and the (distant) torus reflector are constant, we can explain the 
varying fitted reflection strength by the presence of a contribution 
of the flux reflected by the torus, which is constant over the obser- 
vation time scale given its very large size (Section l5.1.2t . In the 
bright state, the contribution of the torus reflection is small, but it 
becomes much stronger in the dim state, explaining the fitted large 
value of R. Given the above assumption, we find that the solid an- 
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gles subtended by the disc and torus (as seen from the X-ray source) 
are - 0.3, 0.2 of In, respectively. These values have been obtained 
with a slab reflection model, which is likely to substantially under- 
estimate the actual solid angle subtended by the torus, which might 
then be close to In. 

The Comptonizing plasma is photon-starved, i.e., the flux in 
seed photons incident on the plasma is ~ 15 times weaker than 
the Comptonized flux (Section 15. 3t . This is consistent with the 
seed photons being supplied by a truncated optically-thick disc sur- 
rounding a hot accretion flow. This geometry is also consistent with 
the relatively small disc reflection. We also find that the flux of the 
disc photons fitted to the X-ray spectra is an order of magnitude be- 
low the UV fluxes actually observed. This is again consistent with 
the above geometry, in which only a small fraction of the emitted 
disc photons is incident on the inner hot regions. All these findings 
rule out a static disc corona geometry. However, they are still com- 
patible with a corona outflowing at a mildly relativistic speed, see 
SectionEOI 

In any case, the disc inner radius cannot be very large given 
that the observed EF E at UV are of the same order as those at hard 
X-rays, implying a rough equipartition between the disc and hot 
flow emission. For a non-rotating black hole, the radius with the 
same integrated dissipation below and above it is ^ \5R„. The trun- 
cation radius can also be constrained by the fitted Comptonization 
model, which assumes the seed photons are disc blackbody. This 
gives a similar value of R m for the inner disc temperature of ~ 5 
eV, which value is consistent with the observations. 

Compared to Seyfert Is, NGC 4151 has a somewhat harder 
X-ray spectrum and less reflection. Its spectrum is also similar to 
black-hole binaries in the hard state, in particular to GX 339^1, 
which has shown a very similar spectral evolution. 

Future missions, e.g., Astro-H and EXIST, with instruments 
with a sensitivity comparable to that of OSSE, can provide better 
data above 200 keV than obtained up to date. Yet, for a better un- 
derstanding of the physics of the central engine of AGNs, we need 
high-quality data up to at least 1 MeV. This would yield strong con- 
straints on the electron distribution in the Comptonization region 
(e.g. the presence of a non-thermal tail) and on the source geome- 
try (e.g. the spatial distribution of the electron temperature and the 
characteristic optical depth. 
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